Spin-charge filtering through a spin-orbit coupled quantum dot controlled via an 

Aharonov-Bohm interferometer 
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We show that a strongly correlated quantum dot embedded in an Aharonov-Bohm interferometer 
can be used to filter both charge and spin at zero voltage bias. The magnitude with which the 
Aharonov-Bohm arm is coupled to the system controls the many-body effects on the quantum dot. 
When the quantum dot is in the Kondo regime the flow of charge through the system can be tuned 
by the phase of the Aharonov-Bohm arm, ipab- Furthermore when a spin-orbit interaction is present 
on a Kondo quantum dot we can control the flow of spin by the spin-orbit phase, ipso- The existence 
of the Kondo peak at the Fermi energy makes it possible to control the flow of both charge and spin 
in the zero voltage bias limit. 

PACS numbers: 73.23.-b, 85.35.-p, 85.75.-d 



I. INTRODUCTION 

The ability to easily control charge and spin trans- 
port is of great importance in nanotcchnology, specifi- 
cally spintronics. The Quantum Dot Aharonov-Bohm 
Interferometer (QD-ABI) (Figure [TJ has been found to 
be a candidate for manipulating electron spins^. In re- 
cent years a number of theories^ have been put forth to 
take advantage of the interference effects in such a geom- 
etry. Sun and Xie^ showed that the spin polarization on 
the QD may be controlled via the voltage bias. In the 
presence of a local coulomb interaction on the QD Hof- 
stetter et alA studied the dependence of the Fano line 
shape on the AB phase. 

In this paper we show that the addition of a spin-orbit 
(SO) interaction on the strongly correlated QD allows 
the QD-ABI to function as a spin-charge filter. The abil- 
ity of this system to act as a filter is a consequence of 
the interference between the continuum (AB arm) and 
the localized state (QD) , which is known as the Fano 
effect^. Aside from the QD-ABI the Fano effect has been 
observed in a single electron transistor—, a quantum wire 
with a side coupled QDL&, and multiwall carbon nan- 
otubes in a crossed geometry^. 

The presence of the AB arm tends to localize electrons 
on the QD, depending on how strongly the arm is coupled 
to the system. Therefore this property will give us an 
effective control over the strength of these correlations 
on the QD. We may exploit this tuning of the many- 
body physics to control the charge-spin transport of the 
QD-ABI system. 

The QD interactions we will study in this paper are an 
on-site Coulomb interaction and the Rashba spin-orbit 
interaction^. The Coulomb interaction in the Kondo 



regime will allow us to filter both charge and spin at 
zero bias. The reason for this is due to the fact that the 
Kondo effect induces a sharp resonance in the QD spec- 
tral function at the Fermi energy. The further addition of 
the SO interaction, which is induced by the application 
of a gate potential, for a single orbital QD, will create a 
spin dependent phase factor— in the AB arm tunneling 



coefficient. 

Often the spin-orbit interaction is considered a cou- 
pling with spin degrees of freedom mediated by interlevel 
transitions in quantum dot systems and therefore, due to 
the significant level spacing and the Coulomb interaction, 
the inter-level SO coupling strengths in QD systems are 
thought to be small. However, as pointed out in Ref. 2 , 
an intra-level phase factor induced by the SO coupling 
may be realized in high g-factor systems^ such as In- 
GaAs quantum dots. Without the interference effect, eg. 
without the AB-arm, such effect can be ignored. How- 
ever, the presence of the AB-arm not only manifests the 
intra-level SO interference effects but also effectively con- 
trols the many-body effects on the QD, hence a strong 
influence on the charge-spin transport. 

The experimental setup of the QD-ABI is consistent 
with that of Kobayashi et al*£. The left and right leads 
are modeled as infinite non-interacting electron reser- 
voirs. They are connected to each other via two arms. 
The top arm is the AB arm which has a complex tun- 
neling coefficient to — \to\e lipAB that is controlled by a 
magnetic field. The sign of tpAB is positive for electrons 
traveling from the left reservoir to the right reservoir. We 
consider the magnetic field to be small enough so that we 
may ignore the Zeeman splitting of the QD energy level. 
The bottom arm contains the embedded QD with real 
tunneling coefficients fx and tn. We choose fj,, tn to 
be real because the L, R states are defined to absorb 
the phase factor. Our calculations are carried out in the 
low-bias, linear response regime. 



II. THEORY 

A. Transmission coefficient T(e) 

In this section of the paper we will derive the gen- 
eral transport functions for the QD-ABI, which are cor- 
rect with the interaction on the QD. The non- interacting 
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FIG. 1: Experimental setup of the quantum dot embed- 
ded AB interferometer. Two infinite electron reservoirs (Left, 
Right) are connected by two arms. The bottom arm contains 
the embedded QD with real tunneling coefficients th and tR. 
The upper arm is a direct connection between the left and 
right leads with a complex tunneling coefficient to = \to\e lIfiAB , 
where ipab is a phase factor controlled by the magnetic field, 
B. The sign of i{>ab is positive for electrons traveling from 
the left to the right, as the arrow indicates. It, Ifl, and lm n g 
are the Left, Right, and Ring currents respectively. 



Hamiltonian of the system is given by 
w = 7~l LR + li d +T~k- 



The Hamiltonian Ti consists of three parts: Hl.r de- 
scribes the reservoirs, H.d the QD, and TL t the tunneling 
between the reservoirs. 
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where c^ fco .(c Q feo-) and d^d^) are the creation (annihi- 
lation) operators with momentum k and spin a of the 
a = (L, R) reservoir and the QD, respectively. In addi- 
tion, O is the volume of the reservoirs which is taken to 
infinity. 

Before presenting the Landauer formula for the cur- 
rent we define the parameters. In the non-interacting 



limit without the AB arm, the line broadening of the QD 
spectral function due to the leads is T = Tl + Tr where 
r a = irNota- In our calculations we take the density of 
states, No, to be a constant. 

The exact current from the L to R reservoir, regard- 
less of the local interaction on the QD, is given by the 
Landauer formuk 
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Il = — J T(e)Af(e)de. 
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Here A/(e) = /k(e) — an d T(e) is the transmis- 

sion function. The transmission function was previously 
reported in Refs — although without derivation. There- 
fore we present the derivation, which makes use of stan- 
dard Keldysh Green function technique s 14 ' 15 , in the Ap- 
pendix. Here we summarize the results. 

The transmission function may first be decoupled into 
two parts, the flow of current through the QD and 
through the AB arm, 



(1) T(e)Af = i QD (e)+i AB (e) 
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i AB (e) = [-t* G< L (e)+t G< R (e)]=2Re[-t* G< L (e)l 



where we have used the relation, GS 



and i AB are the contributions to the current from the QD 
and the AB arm respectively. To simplify our notation 
we define the following parameters 
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where ro = ir 2 N^t^ 2 . Here Tq is the transmission func- 
tion when the QD is disconnected from the left and right 
reservoirs, i.e. tz, = tR. = 0. The current through the QD 
and AB arm is found below. 
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In our analysis we find there exists a ring current even 
at zero bias, 
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Here / = (/l + /fl)/2 is the average Fermi function. This 
current flows in the clockwise direction through the AB 
ring and persists even at zero bias since it is proportional 
to /, not A/. Although the ring current can be of the 
same order of magnitude as the total current, it does not 
contribute to the source-drain current. 

Using Eq. © we arrive at the exact transmission func- 
tion, 



T(e) = T - 2f V^Ro cos(^ B )Re[G« (e)] 
- f HI - T cos 2 (^4 S )) - T ]Im[G« (e)]. 



(14) 



We emphasize that G^ d (e) is the full interacting QD re- 
tarded Green function, and Eq. [14] applies to systems 
with an interacting QD. 



B. Non-interacting limit 

In the non-interacting limit we will address two impor- 
tant points. First, as the coupling to the left and right 
reservoirs is increased through \t \ 7 the electron becomes 
more localized on the QD. Secondly, the noninteracting 
system is not suitable for controlling both charge and 
spin transport. 

The noninteracting QD spectral function is given by 



A(e) = 



t/tt 



where 



S = — — cos(^ AB ). 



(15) 



(16) 



and f was defined in Eq. (TTU)) . We plot A(e) in Figured] 
for different values of |to|. Notice that as the magnitude 
of \to\ is increased, the QD spectral function becomes 
sharper and the center of the peak is shifted. This shift 
in the QD energy, 5 and the reduced line broadening, 
r are easily understood by performing a bonding-anti- 
bonding transformation of the leads in real space. Upon 
doing so the QD is now only coupled to the first site of 
the bonding chain where the local energy of this site is 
shifted by — |io| cos(^ab)- Thus by increasing \ta\, the 
connection of the QD to the bonding chain is reduced, 
and the energy level of the QD is shifted upward. 

The differential conductance is given as G( ( f>) = 
e(dl/d&), and in the zero bias and low-temperature limit 

G(0) = =f-T(0). Therefore at equilibrium we only need 
T(0) to determine the conductance. In Figure [3] we plot 
the transmission amplitude as a function of the AB phase, 
<p AB . When ip. AB = ±f , T(0) = 1. When ip AB = (0,tt), 




FIG. 2: QD spectral function for different magnitudes of \to\ 
at ipAB = 0. As |to| is increased the electron becomes more 
localized on the QD and the distribution of energies shift away 
from the Fermi energy. 



the Fano anti-resonance becomes most prominent and the 
minimum of the peaks are positioned significantly above 
and below the Fermi energy. As a result of this fact we 
are not able to extinguish the charge or spin conductance 
at zero bias. If we were able to shift the minima of these 
anti-resonance peaks to the Fermi energy, then we would 
be able to fully control the transport through this system 
at zero bias by tuning <pab- In the next section of this 
paper we show that we can accomplish this through the 
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FIG. 3: Noninteracting transmission function for different 
values of <pab- When ipAB = (— f, f-) we see a resonance at 
the QD energy level. The ipAB = (— f,?) curves coincide, 
away from these curves the interference effects become more 
pronounced and the transmission function becomes asymmet- 
rical. When ifAB = (0, n) the interference effects become most 
pronounced and a very strong anti-resonance emerges which 
is shifted to the left or right of the QD energy level. In the 
limit |e| — * oo the transmission converges to a finite value, 
due to the AB arm. 



C. Interacting Fano Effect 

When a local many-body interaction is present on the 
QD site, the problem is essentially reduced to calculat- 
ing the full QD retarded Green function, G dd (e). Once 
we know this Green function we may then calculate the 
transmission function using Eq. (|14[) . In this paper we 
take into account the Coulomb interaction on the QD in 
the form of the Anderson interaction, 



'Hint = Uhd^hdi- 
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We perform the diagrammatic calculation in the 
imaginary-time formalism. The QD Green function, at 
imaginary Matsubara frequency iu> n — j ( 2 "+ 1 ) ,r ; j s given 

by G dd (ioJ n ) = [(G% d )~ 1 (iu n )-'E(iu}n)]~' L > where the self 
energy, S(iw„) is calculated to second-order in U (Fig- 
ure |4|) . The second-order perturbation theory has been 
studied extensivel y 17 ' 18 i 19 i 20 i 21 and shown to be a very 
good approximation in the particle-hole symmetric limit 
up to values of U/T ~ 6 (weak-coupling regime) . This 
weak-coupling approximation has also been used within 
the framework of dynamical mean field theory^, and 
produced agreement with nonperturbative methods. In 
our model the particle-hole symmetry will be broken for 
5^0, but the calculation of the zero bias conductance is 
in excellent agreement with the nonperturbative numer- 
ical renormalization group^NRG) as will be shown. 
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FIG. 4: Self energy expanded to second order in U. 

The first-order diagram becomes U (rid) = ^ hi the 
half- filled limit. Absorbing the first-order diagram into 
the non-interacting Matsubara Green function we have 
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and the second order self energy becomes 
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Now we are in position to solve this problem numer- 
ically. First we calculate the self-energy in Matsubara 
frequency. Then in order to calculate T(e) we must ana- 
lytically continue the Green function to its retarded form 
in real frequency space, Gdd{iu n — * e + irf). For the 
numerical analytical continuation we used the A-Point 
Pade approximant method 2 ^. 

Let us first look at the spectral function, A(e) = 
— —lm[Gdd(c)]- For the resonant case (ipAB — §) the 




FIG. 5: Interacting QD spectral function for different mag- 
nitudes of |to| when 95 ab = ■§, U = 0.6 and /3 = 160. We see 
that as \to\ is increased the electron becomes more localized 
on the QD and as a result the correlations due to the Coulomb 
interaction become more pronounced. 



spectral function is plotted for different values of |to| hi 
Figure O We see that as we increase the coupling |t | of 
the L, R states to the AB arm, the effective many-body 
interaction U/T is strongly enhanced. From our analysis 
of the noninteracting system this is exactly what we ex- 
pected to happen. When \t \ =0 we start in the weakly 
interacting valence fluctuating regime, but as we increase 
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FIG. 6: QD spectral function and transmission amplitude 
for U = 0.6 and ft = 160. The if>AB = f, -rr curves are iden- 
tical, and at these values the spectral function is symmetric 
and the transmission amplitude displays resonant behavior. 
On the other hand when ipAB — 0, n the spectral function 
is asymmetric and the transmission amplitude has an anti- 
resonance. 

The transmission functions is given in Figure [5] As 
in the noninteracting case we have resonance (anti- 
resonance) phenomena at (fAB — — -§,-f(0, 7r) respec- 
tively. The anti-resonance peaks are known as the Fano- 
Kondo anti-resonance and have been observed experi- 



5 





i 1 i yv i 1 




\\ // M Hol/tL/R = 1 




\\^y / °-°ltol/tivR=1.40 e 




\ / — lt l/t L /R = 1.90 




\ J ^lt l/t L/R = 2.40 




\<*y n= 



<Pab 

FIG. 7: Zero bias conductance as a function of tfAB for 
th = tR = 0.15, U = 0.5 and /3 = 160. The conductance 
approaches unity when tfiAB = 2 " 2 +1 7r. Conversely the con- 
ductance approaches a minimum when ifAB = mr, with inte- 
ger n. 



mentally in a quantum wire with a side coupled QD 7 - 8 . 
For ifAB = 0, 7r the anti-resonance peaks are antisym- 
metric of one another about the Fermi energy. 

Now let us examine the zero bias conductance. Our 
results are given in Figure [7] as a function of fAB for 
different values of the ratio \to\/tL/R- The NRG calcu- 
lation of RefjS' examined the zero bias conductance as a 
function of fAB for different values of the gate potential, 
e d . A careful straightforward identification of the model 
parameters shows that our results agree excellently with 
the NRG 25 , which justifies our self-energy approxima- 
tion. Since G(0) is only dependent upon the value of the 
Green function at the Fermi energy, this implies that the 
second-order self energy approximation at least produces 
reliable results near the Fermi energy for S =/= 0. 

From Figure [7J we see that G(0) oscillates as a func- 
tion of <pab and the magnitude of these oscillations is 
strongly dependant on the ratio, \to\/t L / R . We find that 
\to\/tL/R = 1 is a special case which maximizes the mag- 
nitude of the AB oscillations. These oscillations are 
a consequence of the cos 2 (fAs) term in Eq. (fl"4|) and 
the sharp feature near the Fermi energy of the spec- 
tral function (Figure The resonance (G(0) = max) 
[anti-resonance(G(0) = min)] peaks occur at fAB = 

2 7r[n7r], with integer n. Thus by tuning <pab we may 
filter the charge through the system at zero voltage bias. 
The reason why we have so much control over the con- 
ductance at zero bias is due to the Kondo effect. The 
Kondo effect induces a sharp peak near the Fermi energy 
in the spectral function, insensitive to fAB ■ This feature 
is due to the many body effect and is not present in the 
non-interacting case (Figure [3}. 



D. Spin Transport 

For a single orbital QD the addition of the spin-orbit 
interaction induces a spin dependent phase in tR — , i.e. 



tR — > tRe~ lc "- pso , where a =+(-) for spin up(down). 
Unlike fAB which arises from the orbital motion, fso 
depends on the spin. To simplify our Hamiltonian we 
make the unitary transformation 



so that 



-Rk ^ c c Rk 



to 



(20) 



5 lsvso t (21) 
Now we define a spin dependent AB tunneling coefficient 

\ tQ \ e i(<PAB+*<PS0) _ (22) 



tr, 



As a result of the spin dependence in too-, the Green's 
functions also gain a spin dependence and the spin de- 
pendant second order self energy is given by 
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Now let us examine the spectral and transmission func- 
tions given in Figure [5J As in the spin independent case 
we see resonance and antiresonance behavior in the trans- 
mission function. More importantly, due to the SO in- 
teraction we now have two phase factors, tpAB and tpso, 
which along with the Kondo resonance near the Fermi 
energy, we may use to filter the spin up and spin down 
electrons independently. 

Let us take a closer look at Figure [H We see that when 
fAB = fso = f > -A-i has a spectral structure similar to 
that of a Kondo dot while A± develops a slight asymme- 
try from the spin up case. Looking at the transmission 
function, Tf is similar in shape to its spectral function; 
and Tf shows strong antiresonance behavior due to the 
AB ring. This mechanism gives us a strong control on 
the spin-transport. At the Fermi energy Tj(0) = 1 and 
Tj. = 0. When ipso = f both T T and T± show strong 
interference effects. 

To look at the zero bias spin conductance we turn to 
Figure [9] where we show the SO oscillations of the zero 
bias conductance. The behavior of the spin up and spin 
down conductance is a result of the cos 2 (tpAB + Vso) 
term for spin up and the cos 2 (<pab — <pso) term for 
spin down in Eq. (|14[) . In the figure we see that when 
fAB = f the spin polarized conductance [77 = (Gf (0) — 
G|(0))/(G|(0) + Gx(0))] is maximized. Conversely when 
if so = f the spin polarization is suppressed. 

The most important results of this paper are shown 
in the spin polarized conductance in Figure 1101 Here 
we plot 77 as a function of ipso f° r different values of 
fAB- The maximum in the spin up/down conductance 
occurs approximately when fAB ± fso = 2n 2 t " 1 7'"; while 
the minimum occurs at fAB ± fso — Tnm. This can 
be seen from the cos 2 (y) factor in Eq. (fl"4|) , since at the 
Kondo resonance (e « 0) only To and the term propor- 
tional to lm[G dd (e)] become relevant. Therefore the nec- 
essary conditions on fso and fAB for the spin polarized 



conductance to be maximized and minimized are 
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n — m + i 



n + m + 



n + m + 1 



(24) 
(25) 



Here again we emphasize the crucial role which the 
many-body interactions play in this system. The 
Kondo effect serves to simultaneously pin the reso- 
nance/antiresonance peaks of the spin dependent con- 
ductance at the chemical potential, as seen in Figure [5] 
As a result, when the SO interaction is turned on, we gain 
complete control over the spin polarization (Figure [TTJ]) . 
Furthermore, even if there only exists a small SO inter- 
action in the QD, this device serves to enhance those 
inherent SO effects. In contrast to this, in the nonin- 
teracting device (Figure [3]) the resonance/antiresonancc 
peaks occur at significantly different energies, making it 
impractical to use as a spin/charge filter. Therefore we 
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FIG. 8: Spin dependant spectral and transmission func- 
tions for ipso = f,fi U = 0.6 and (3 = 160. For ipso = 
? Af(Ai) is symmetric(asymmetric), while T|(Tj) displays 
resonance(anti-resonance) behavior. In addition, T\(0) ~ I 
and Tj.(0) ~ 0. For ipso = § both Af and Ai axe asymmet- 
rical and A i lies on the curve of A \ (ipso = j)- and Tj, 
both display anti-resonant phenomena. 



III. CONCLUSION 

In this paper we have shown that the QD-ABI may be 
used to localize and delocalize electrons on the QD. As a 
result, when the QD is strongly correlated we can effec- 
tively control the strength of these correlations through 
the magnitude of to- With larger to, electrons get more 
localized on the QD. In the case of a coulomb interaction 
on the QD, the Kondo effect induces a sharp peak in the 
spectral function near the Fermi energy. We emphasize 
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FIG. 9: SO oscillations for ip AB = £, f , U = 0.6, and (5 = 
160. The minima and maxima of the SO oscillations approach 
and I respectively. 
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FIG. 10: Spin-polarized conductance for different AB phases 
when U = 0.6 and /3 = 160. Here the spin polarization is 
defined as, 77= [Gf (0) — Gj.(0)]/[G T (0) + Gj.(0)]. We see that 
the spin polarization may be fully controlled by tuning f>AB 
and (pso- The ipAB = 0, ^ curves are identical. 



that the exact relation for the transmission, Eq. (|14|) . 
shows that the resonance/anti-resonance is driven by the 
AB phase factor once the spectral function develops a 
sharp feature near the Fermi energy, independent of its 
many-body character. 

The most important result of this paper is realized 
when a SO interaction is present on the QD in addi- 
tion to the coulomb interaction. In this case the Kondo 
peak and the spin dependent cos 2 ^) term induce the 
resonance/anti-resonance in the zero bias spin dependant 
conductance. Due to this property we can vary the Gf (0) 
and Gj(0) between and e 2 /h. Further, the SO phase 
gives us another degree of freedom, in addition to the AB 
phase, so that the spin up and spin down conductances 
may be varied independently. Thus by tuning <pab and 
'Pso the spin polarization may be fully controlled in the 
QD-ABI. 
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APPENDIX A: DERIVATION OF T(e) 



and define the following QD-excluded Green functions 



In this appendix we present the full derivation of the 
transmission function. The transmission is given in terms 
of Keldysh Green functions, 



T(e)A/ = i QD (e)+i AB (e) 



(0 



-t L [G, 



(Al) 
(A2) 



i AB (e) = [-t* G< L (e)+t G< R (e)] , (A3) 

where G dL and G RL are the Fourier transforms of the fol- 
lowing time dependent nonequilibrium Green functions 
(NEGFs), 



G&(t) 



kk' 



and G< d 



(G<)* and G< 



LR 



<g< l y 



(A4) 
(A5) 
In calcu- 



lating these Green functions it becomes convenient to 
disconnect the QD from the reservoirs, i.e. ti, = tn = 0, 



9ll 



9rl 



9ll 



9rr 
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9rr 
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(fL + r f R ) 



(fR + roh) 

-2tt 2 N$ 
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to Ul 



(A6) 
(A7) 
(A8) 
(A9) 
Ir) • (A10) 



We express G dL and G RL in terms of the fully interacting 
retarded and advance QD Green functions and the above 
QD-excluded noninteracting NEGFs. 

To simplify our notation we define F^ d , F£ d , and F< d 
which are the retarded, advanced, and lesser Green func- 
tions which describe the transport from the QD to the a 
lead in terms of the QD-excluded Green functions. 



F r Ld = (FSlT = h (t L g LL + t R g LR ) r 



F< 

r Ld 



r Rd 



1 



l + r 
1 



F r Rd = (FS R T = t R (t R g RR + t L g RL ) r = —— 1- % T R + v^TTT^e 



l + r 



~(F< R )* = t R (t R g RR + t L g RL )< = — — \iT R (f R + r f L ) + v^IV^ (f R - f L j 



(All) 
(A12) 
(A13) 
(A14) 



r 



Therefore G RL and G dL may be written in terms of F's 
and the full QD Green function as 



t Lt rG R l — F Rd G dd F dL 



~t L G dL — G dd F dL . 



Using the Keldysh Green function relation: 



,2(, 



(AB)< = A < B a + A r B < 
(ABC)< = A<B a C a + A r B<C a + 
A r B r C<, 

the lesser Green functions become 



Grl 



-IlG^l 



Z7< y^a rpa I 
^Rd^dd^Ld + 

t Rd (j dd t Ld + r Rd Lr dd- t Ld 

dd^dL + ^dd? dL- 



(A15) 
(A16) 



(A17) 
(A18) 



(A19) 
(A20) 



Making use of the nonequilibrium steady state condition, 
\l + i R = 0, we may construct G dd in terms of G r dd and 
G dd as follows. The ensemble averaged currents are given 
by 



- -t Q (G da - Gl 
= (F da ~ F^ d )G d 



U dd^da 



F< d G a dd (A21) 



where a = (L, R). Therefore by invoking the steady state 
condition the QD lesser Green function becomes 



Gdd 



(F ( < 



dL 



F d < R )G 



dd 



( F Ld + F Rd )G a di 



<ld 



(F£d F dL ) + (F Rd F dR ) 



(A22) 



Inserting G dd into Eq. (|A21|) with a = L we arrive at 
Iqdi Eq. (fTTj) . The current through the AB arm is given 
by 



8 



= {-toFm F dL + t O F Ld F ^d^dd + (~ t: F Rd F dL + ^Ld^R^dd 



{-ta F Rd F dL + to F Ld F dR)G. 



dd • 



where all of the terms have been solved for. Thus after 
doing some algebra, i AB may be expressed as Eq. (fT^j). 
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